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The  delamination  of  orthotropic  laminates  containing  finite-length  cracks  and  subject  to  thermal  gradients  is  examined. 

The  exact  limiting  case  solutions  for  infinitesimal-  and  infinite -length  cracks  are  known,  and  are  equal  to  each  other  when 
the  crack  length  is  approximately  equal  to  the  plate  thickness.  However,  in  the  transition  region  of  crack  length  from  about  1 
to  5  times  the  plate  thickness,  both  limit  solutions  overestimate  the  energy  release  by  20-100%.  Hence,  an  analysis  was 
developed  to  better  predict  the  energy  release  rate  for  such  finite-length  cracks.  The  model  is  a  modification  of  the  infinite- 
crack  analysis  of  Hutchinson  and  Lu  (1995,  ASME  J.  Eng.  Mat.  Tech.,  117  (4)  pp.  386-390)  and  provides  a  closed  form 
expression  for  the  elastic  energy  release  rate  in  a  plane-strain  orthotropic  flat  plate  that  agrees  well  with  numerical  values  for 
cracks  of  length  approximately  half  of  the  plate  thickness  and  larger.  The  analytic  result  is  shown  to  agree  well  with  finite 
element  results  over  a  wide  range  of  crack  lengths,  depths  and  interface  conductivity,  both  for  isotropic  and  orthotropic 
materials. 


Introduction 

Advanced  high  temperature  materials  such  as  ceramic  composites  are  under  active  investigation  for 
use  in  hot  shell  structures  such  as  combustion  liners.  These  materials  are  typically  characterized  by  low 
thermal  conductivity  and  a  laminar  structure  with  relatively  low  interlaminar  strength.  Consequently, 
the  hot  structures  will  be  subject  to  large  thermal  gradients  and  thermal  stress  induced  delamination 
initiated  at  internal  interlaminar  flaws  are  of  significant  concern.  The  thermal  delamination  problem  has 
been  studied  extensively  in  the  literature.  The  stress  concentration  due  to  an  infinitesimal  crack  in  an 
elastic  body  subject  to  a  thermal  gradient  is  given  by  Sih  [1],  The  complementary  asymptotic  linear 
elastic  fracture  solution  for  an  infinitely  long  crack  in  an  isotropic  flat  plate  is  given  by  Hutchinson  and 
Lu  [2]  (H-L).  These  basic  solutions  are  followed  in  the  literature  by  analysis  of  more  complex  systems. 
Sorensen,  Sarraute,  Jorgensen  and  Horsewell  [3]  analyzed  the  long-crack  problem  for  a  laminate  with 
dissimilar  elastic  laminae,  and  show  that  the  interfacial  energy  release  rate  may  be  mitigated  by 
selective  layering  of  materials.  Gu  and  Asaro  [4]  considered  delamination  fracture  in  continuously 
graded  materials.  The  elastic-plastic  problem  was  considered  by  Aoki,  S.,  Kishimoto,  K.,  and  Sakata 
[5], 

The  present  work  focuses  on  the  homogeneous  orthotropic  plate  problem.  Numerical  investigation 
reveals  that,  with  increasing  crack  length,  the  transition  from  the  Sih  infinitesimal  crack  result  to  the  H- 

L  infinitely  long  crack,  or  steady-state,  solution  is  less  abrupt  than  estimated  by  Hutchinson  and  Lu. 
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Here,  we  present  a  modification  to  the  H-L  result  that  provides  a  substantial  improvement  for  finite 
length  cracks  (i.e.  on  the  order  of  the  plate  thickness),  while  retaining  the  exact  limiting  result  for 
arbitrarily  long  cracks. 

The  numerical  work  that  motivated  the  present  analysis  was  a  straightforward  2-d  plane  strain 
ABAQUS1  finite  element  model  of  a  central  crack  in  a  flat  plate  subject  to  a  prescribed  temperature 
differential  between  the  surfaces  as  shown  in  Figure  1.  A  focused  mesh  with  second  order  quarter  point 
elements  at  the  crack  tip,  and  reduced  integration  elements  throughout  was  employed[6].  Energy 
release  rates  were  calculated  using  the  thermoelastic  J-integral  formulation  built  into  the  commercial 
code  [7,8].  For  verification,  the  calculated  energy  release  rates  were  compared  with  the  exact  long- 
crack  H-L  solution  and  were  found  to  be  25%  below  the  asymptotic  prediction  for  cracks  as  long  as  3 
times  the  plate  thickness,  with  convergence  to  the  asymptotic  result  confirmed  only  by  modeling  cracks 
20-30  times  the  plate  thickness.  Thus  it  seemed  appropriate  to  develop  an  analytic  model  to  better 
predict  the  fracture  behavior  for  such  intermediate  length  cracks. 

The  H-L  model  is  based  on  a  long-crack  isothermal  fracture  model  of  Hutchinson  and  Suo  [9]  along 
with  several  infinite  crack  approximations  for  the  temperature  distribution  and  resulting  deformation 
due  to  a  thermal  gradient.  Numerical  simulation  of  the  isothermal  model  shows  that  it  performs 
remarkably  well  for  cracks  as  short  as  twice  the  plate  thickness.  There  are  two  key  asymptotic 
kinematic  assumptions  in  the  thermoelastic  analysis.  1)  It  is  assumed  that  the  upper  and  lower  sections 
of  the  plate  should  have  equal  curvatures.  2)  The  in-plane  strain  parallel  to  the  crack  faces  at  the  crack 
center  is  assumed  to  be  equal  on  the  opposed  crack  faces.  By  inspection  of  the  numerical  simulation, 
the  equal  curvature  assumption  proves  to  be  quite  good  for  relatively  short  cracks.  The  equal  strain 
condition  is,  however,  approached  rather  slowly  with  increasing  crack  length,  and  here  an  improvement 
is  suggested. 

Basic  analysis  summary 

The  essential  calculation  is  the  same  here  as  in  H-L,  and  so  here  we  will  highlight  the  key  points  as 
briefly  as  possible.  The  thermal  load  results  in  a  set  of  resultant  forces,  P,  and  moments  Mx ,  M2  per  unit 
out-of-plane  thickness  at  the  center  of  the  crack  as  shown  in  Figure  2.  The  stress  intensity  is  then 
calculated  in  terms  of  those  end  loads  using  the  isothermal  fracture  formulation  set  forth  in  [9]. 
Referring  to  Figure  2,  equilibrium  requires  Mx  +M2  =  PH/ 2  .  For  a  long  crack  the  upper  and  lower 
plate  sections  are  modeled  as  beams  and  the  curvatures  are  asserted  to  be  equal  far  from  the  crack  tip, 
leading  to  the  condition  Mx/Hx  =  M2 / H\  ,  thus 


1  Abaqus,  Inc.,  Pawtucket  R.I. 
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(1) 


M1  = 


rfPH 

2(1  w) 


,  M,= 


PH 

WW) 


where  77  =  Hx  /h2  ■  The  analysis  is  restricted  to  0  <  77  <  1 ,  ie.  if  the  crack  is  off  center  it  must  be  nearer 
to  the  hot  side  of  the  plate  to  ensure  an  open  mode  crack  [2],  Equation  (1)  is  as-given  in  [2],  however 
here  we  retain  the  resultant  force  P  as  an  unknown. 


The  orthotropic  thermo-elastic  constitutive  relation  is  of  the  form, 

=  ^/EX~  VyX°yy  lEy~  Vzz  I E : z  +  « AT 

sxz=(yxJ  2pxz 


where  sVj ,  cr .  are  the  components  of  strain  and  stress  and  expressions  for  the  remaining  strain 
components  are  obtained  by  suitable  change  of  subscripts.  Ei,juij,vij  are  the  orthotropic  elastic 
constants  which  satisfy  the  symmetry  requirement,  vijEj  =  vjiEi ,  a;  are  orthotropic  thermal  expansion 
coefficients  and  AT  is  a  change  in  temperature. 


Using  Eq.  (1),  and  assuming  plane  strain  {e  =  syz  =  s  =  0)  the  stress  intensity  factors  are  found 
directly  using  the  analysis  given  in  [9].  In  the  centered  crack  (7  =  1)  case  the  crack  tip  field  is  pure 
mode  II,  with  stress  intensity  factor, 


K 


7-1 


\II 


1/2 


(3) 


where 


A  = 


1  +  P 

2/l1/2 


\!/2 


(4) 


X  and  p  are  the  orthotropic  material  constants  defined  in[9]  and  in  Lu,  Xia  and  Hutchinson  [10], 


X  =  EjE 


P 


=  T1/2 


V  +V  V 

xz  xy  yz 


2/U  1  j 


(5) 

(6) 


where  E  =  Ex/(l-vxyvJ  and  E_  =  Ez/(\-vzyvyz)  are  the  effective  plane  strain  moduli  for  in-plane 

uniaxial  stressing  in  each  of  the  principle  directions.  See  also  Krenk  [11],  using  p  =  k  and  2  =  S~4 .  Note 
that  X,p,  and  A  are  all  unity  for  an  isotropic  material. 


In  the  general  case  for  0  <  tj  <  1  the  crack  is  mixed  mode  with  the  energy  release  rate  and  stress  intensity 
factors, 
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(7) 


x„=/c;r'r 


k,=X!\geI\-kiY 

Where  T  is  the  ratio  of  mode  II  stress  intensity  to  the  rj  =  1  stress  intensity, 


r  = 


sin  co  +fe+iy  cos  <»  +  y 


(7a) 


2^  J  l2V2[7  4(l +  73)  V2F 

and  U,  V,  w,  y  are  the  dimensionless  functions  of  // , 

i/f/= (1+7X1 +3^(1+^)) 
i/r  =  i2(i+^3) 

sin  ;*  =  6t/2  (l  +  jj\JUv  1  ' 

<y  s  3°(l  - ;/)+  sin-1  2 /  V7 

Kinematic  constraint 


In  place  of  the  H-L  condition  of  strain  equality  at  the  crack  center,  the  kinematic  condition  we  will 
impose  is  that  the  deformed  lengths  of  the  upper  and  lower  crack  faces  must  be  equal.  This  leads  to  the 
condition, 

](s[+J-s^)dr  =  0  (9) 
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where  sxx  is  the  strain  parallel  to  the  crack  face,  r  is  the  coordinate  along  the  crack  measured  from  the 
tip,  and  superscripts  (+)  and  (-)  denote  the  upper  and  lower  faces  respectively.  The  orthotropic  plane 
strain  constitutive  relation  on  the  traction- free  faces  is, 


,(±f 


r(±) 


+aT 
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TrT^) 


(10) 


where  cr^  is  the  stress  parallel  to  the  crack  face,  a  =  ax  +  v yx a y  ,  is  the  plane  strain  thermal  expansion 
coefficient,  and  T  is  the  crack  face  temperature  The  (±)  notation  is  used  to  indicate  that  this  is  a  pair  of 
equations  separately  relating  the  quantities  on  each  face.  Combining  Eq.  (9)  and  Eq.  (10)  yields, 

\^x) -<j{h})dr  +  Eaa<f>{l i -T2)=  o  (11) 

0 

where  T] ,  T2  are  the  temperatures  at  the  plate  surfaces  and  we  define  the  average  temperature  jump  <j> , 

=  (r(+)-r(-})rfr  (12) 

1 1  —  1  o  Cl  ~ 

1  7  0 
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For  comparison  with  the  H-L  formulation,  for  a  sufficiently  long  crack,  far  from  the  crack  tip,  the  heat 
conduction  through  the  plate  is  ID  and  the  jump  in  face  temperatures  away  from  the  crack  tip  will  be 
determined  by  the  conduction  across  the  crack  as, 


rp(+)  _  f(- ) 


T  -T 
1\  12 


1 


x=0 


1+5, 


(13) 


where  the  Biot  number,  Bc ,  characterizes  the  heat  flow  across  the  open  crack,  with  Bc=  0  for  a 
perfectly  insulating  crack,  and  Bc  — >  oo  for  perfect  conduction[2].  Near  to  the  tip,  heat  flux  around  the 
tip  reduces  the  temperature  jump  across  the  faces.  However  in  the  limit  a»  H  the  influence  of  the 
crack  tip  field  becomes  vanishingly  small,  so  that  ^  » (l  +  5C  )_*  •  For  any  finite  length  crack 

0<^<(l  +  5c)-\ 


By  inspection  of  numerical  results,  for  cracks  that  are  approximately  equal  to  the  plate  thickness  and 
larger,  the  stress  along  the  crack  faces  is  well  approximated  by  matching  the  asymptotic  crack  tip  (Paris 
and  Sih,  [12])  and  far- field  (beam  theory)  solutions, 


cr*Y  =  < 


r(±) 


r*(±)  <  r  <  a 


-s(±)Ku  A| 


Y/2 


\nrj 


0  <r  <r 


i±) 


(14) 


where  cr^  is  the  value  of  the  tangential  stress  at  the  plate  center,  and  the  sign  indicator 
=  UY-1  =  -1  ,  is  introduced  because  the  stress  is  compressive  on  the  upper  face  and  tensile  on  the 
lower.  The  transition  radius,  r*^>  is  the  distance  from  the  crack  tip  (generally  different  on  each  face) 
where  the  two  limiting  forms  are  equal,  which  is  simply, 


i±)=l 

71 


A K 


II 


r(±) 


; 


(15) 


Figure  3  shows  the  piecewise  stress  approximation  compared  with  numerical  results.  It  should  be  noted 
that  the  piece-wise  function  is  somewhat  lower  than  the  actual  distribution  near  r*^\  which  ultimately 
results  in  an  over  (i.e.  conservative)  estimate  of  the  energy  release  rate.  The  crack  face  stress  far  from 
the  crack  tip,  ajf1 ,  is  found  directly  from  beam  theory  (as  in  H-L), 

cr(1)  =  _i(,)^g(±)(7)  (16) 

M 


where, 

g[±\v) 


(!  +  vf  f,  .  „(±)  (1  -  ^7)(l  -  4/7  +  ?72 

(l  +  7)3  J 


(17) 


so  that, 
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r  ^  V 


4  n 


KSt±]J 


(18) 


recalling  that  T  is  the  stress  intensity  ratio  given  by  Eq  (7a).  Clearly  the  present  analysis  must  be 
restricted  to  cases  where  a  >  r*{±} .  However  as  r*(±)  represents  the  size  of  the  crack  tip  dominated 
region  we  will  conservatively  also  restrict  application  to  Hx  >  2r  (±\  which  is  ensured  for  the  isotropic 
case  if  0.2  <  77  <  1 .  The  maximum  value  of  r*^  in  that  range  is  ~  H/ 10  at  t]  =  0.5 . 


The  crack  length  condition,  Eq.  (11),  can  now  be  integrated  and  solved  for  P, 


1  + 


HA  T2 

Am  g^g^ 


(19) 


In  the  large  a  limit,  and  using  (j)  =  (l  +  Bc )  1 ,  we  recover  the  H-L  result  exactly, 


P  =  - 


(l  +  7)5(l  +  JBc) 


(20) 


Finally,  the  stress  intensity  and  energy  release  rate  are  found  by  substituting  Eq  (19)  into  Eq  (7),  e.g. 


G 

EH[a(T,-T2jf 


V 


(1 +  73V2 


(l  +  7)52 


1  + 


HA 


Am  gK  g 


(+)„(-) 


(21) 


It  is  notable  that  the  quantity  T2/ gi+}gi~>  is  approximately  1  for  0.5  <  //  <  1 .  The  large  a  limit  of  (19)  of 
course  agrees  with  H-L,  hence  we  will  use  the  a  — »  00, 77  =  1  ,  perfectly  insulating  (</>  — >  l) ,  case  as  a 
baseline  for  comparison  of  solutions, 

GHL=EH[a{Tl-T2)]2/ 32  (22) 

For  reference  the  small-crack  isotropic  limit  of  Sih  [1]  is, 

G0  =  27tGHL  j  (23) 


Temperature  distribution 

As  will  be  shown,  Eq  (19)  represents  an  improvement  over  Eq  (20)  even  if  the  crack  face 
temperature  distribution  is  approximated  by  taking  </>  =  (l  +  Bc )  1 .  However,  further  improvement  is 
realized  by  calculating  </>  from  Eq.  (12)  using  the  true  temperature  distribution.  Unfortunately,  although 
the  asymptotic  temperature  distributions  near  the  tip  and  far  field  are  known  (Florence  and  Goodier, 
[13]),  assuming  a  piecewise  temperature  distribution  similar  to  Eq.  (14)  does  not  prove  sufficiently 
accurate.  Therefore  we  will  rely  on  numerically  generated  results  to  derive  an  functional  dependence 
of  $  on  the  plate  geometry.  The  thermal  finite  element  analysis  was  performed  using  a  thermal 
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contact  relation  over  the  crack  surface  such  that  the  heat  flux  across  the  crack  is  q  =  Bck[r(+)  , 

where  k  is  the  (isotropic)  thermal  conductivity  of  the  plate.  Note  that  the  heat  transfer  across  a  physical 
crack  is,  in  general,  strongly  dependant  on  the  crack  opening,  however  here  we  have  assumed  that  local 
variation  of  the  crack  opening  over  the  crack  length  may  be  neglected  and  that  Bc  is  an  effective  value 
which  accounts  for  the  average  crack  opening.  The  consequences  of  this  assumption  will  be  discussed 
in  a  subsequent  section. 


A  high  quality  general  fit  to  the  numerical  results  is  found  by  arguing  as  follows.  For  sufficiently  long 
cracks,  increasing  the  crack  length  does  not  significantly  affect  the  non-uniform  field  near  the  tip, 
further,  for  an  infinite  crack  we  require  </>  =  (l  +  Bc  f 1 ,  therefore  we  fit  our  numerical  simulation  results  to 
the  form, 


<*  = 


1 


(24) 


where  C  is  a  function  ofi?e  and  77.  C  is  determined  from  numerical  solution  of  the  heat  transfer 
problem  as  the  large  a  limit  of  the  quantity  (l  -  (l  +  Bc  )tp\a/H ) ,  which  converges  to  a  constant  quickly 
with  increasing  crack  length  ((a/H)  =  2-3  is  sufficient).  We  find  that  for  cracks  longer  than  a  =  H , 


1 

i  +  B~c 


M 

(l +  Bc  )2/3 


\ 

H_ 

a 


(25) 


where  the  dependence  on  Bc  is  empirical.  The  value  of  £  is  only  weakly  dependant  on  77  over  the  range 
0.5  <  77  <  1 ,  having  values  of  £=;0.2112  at  77  =  1  and  £  =  0.2288  at  77  =  1/2.  As  Figure  4  shows,  this  is  an 
excellent  fit  to  numerically  derived  values  for  all  Bc  and  for  2  a>  H  .  Equation  (25)  is  clearly  not  valid 
for  a/H  <  £(l  +  Bc  )“2/3 ,  but  this  restriction  is  in  line  with  other  aspects  of  the  analysis. 

Results  and  Discussion 

Figure  5  shows  the  energy  release  rate  calculated  using  Eq  (21),  with  <j>  from  Eq  (25)  as  well  as  the 
approximation,  ^  =  (l  +  5c)  1 ,  for  comparison.  The  full  calculation  shows  good  (and  as  expected 
slightly  conservative)  agreement  with  the  finite  element  results.  Ignoring  the  near-tip  temperature 
variation  results  in  a  substantially  higher  energy  release  rate  prediction,  although  it  still  provides  some 
improvement  over  the  asymptotic  result. 

Figure  6  shows  the  effectiveness  of  the  analysis  for  variation  of  both  r\  and  Bc.  The  agreement  with 
numerical  results  is  best  for  77= 1  and  somewhat  decreases,  becoming  more  conservative,  for  off-center 
cracks.  As  noted  previously  this  analysis  is  not  intended  for  significantly  smaller  values  of  rt  (as  for 
example  subsurface  cracks),  where  the  difference  undoubtedly  increases. 
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It  is  worthwhile  to  refer  to  Figure  6  in  discussing  the  effect  of  the  constant  Bc  assumption.  The  analysis 
shows  a  slightly  increasing  energy  release  rate  for  decreasing  7/  at  a  fixed  value  of  Bc .  However, 
realistically,  the  heat  transfer  coefficient  itself  decreases  substantially  (potentially  by  orders  of 
magnitude)  as  the  crack  opens  and  the  mechanism  of  transfer  transitions  from  contact  conduction  to  gas 
transport  and  radiation  (H-L).  As  t)  is  decreased,  the  stress  intensity  becomes  increasingly  mode  I, 
therefore  the  crack  opens  further.  Consequently  if  the  dependence  of  heat  conduction  on  the  crack 
opening  is  accounted  for  then  the  dependence  of  G  on  r]  may  be  far  more  pronounced  than  suggested  by 
Figure  6.  This  aspect  of  the  problem  is  discussed  in  detail  in  H-L  and  is  not  qualitatively  affected  by  the 
present  analysis. 

Closure 


In  closing  we  note  that  although  the  examples  given  have  been  for  isotropic  materials  the 
formulation  holds  quite  well  for  orthotropic  materials  as  well.  Eq.  (19)  predicts  that  the  energy  release 
rate  will  be  reduced  for  large  values  of  the  parameter  A ,  while  the  H-L  result  is  independent  of  any 
orthotropy  (the  orthotropy  affects  only  the  mode  mixity).  Large  A  corresponds  to  small  values  of 
through-thickness  moduli,  Ez  and  juxz ,  which  is  notably  typical  of  many  composite  systems.  As  an 
example  we  consider  in-plane  isotropy,  and  fix  vxy  =  v  =0.1  ,  vxz  =  v  =0.15  and  A  =  1/2 .  Figure  7. 
shows  the  theory  along  with  finite  element  results  for  rj  =  1  and  rj  =  1/2  as  the  in-plane  shear  modulus 
is  varied.  The  model  correctly  predicts  the  reduction  in  release  rate  for  low  shear,  however  the  effect  is 
significant  in  this  example  only  for  very  small  values. 
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Nomenclature 


a 

H,HvH2 

V 

x,y,z 

r 

T 

T  T 

P,  mx,m2 

£ij^v 

Ei’Pij’Vij 

<h  - 

E^Ez 

a 

K„Kn 


crack  half-length 

total,  upper,  and  lower  thickness 

thickness  ratio 

orthogonal  coordinates 

coordinate  from  tip  along  face 

temperature. 

temperatures,  plate  surfaces 
resultant  force  and  moments 
components  of  strain  and  stress 
elastic  constants 
thermal  expansion  coefficients 
effective  moduli 
effective  thermal  expansion 
stress  intensity  factors 
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Kn=' 

iV// 

G 

G0,Ghl 
K  pA 
0 

Bc 

<T0 

(±) 

S'-  ' 

* 

r 

k 

q 


stress  intensity  factor  for  77  =  1 
energy  release  rate, 
asymptotic  limits  of  G. 
orthotropic  material  constants 
average  temperature  jump 
Biot  number 
far  field  crack  face  stress 
sign  indicator 
transition  radius 
thermal  conductivity 
heat  flux 


dimensionless  functions  of  7/ 


u,v,w,r 

g 

r 

£ 

superscripts 
(+)  (-) 

(±) 


stress  intensity  decomposition 
Eq  (17),  g(jl  =  l)  =  1 
ratio  K „(}])/  Kp  ,  t(tj  =  l)  =  1 
fit  parameter 

upper  and  lower  crack  face 

quantity  in  equation  applied  separately  to  each  face. 
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Figures 


Figure  1 .  Geometry  of  an  infinite  flat  plate  with  a  finite  length  central  crack.  Configuration  shown  is  a  =  H ,  Hxj  H2  =  1  /  2  . 
Contour  shading  shows  the  numerically  calculated  temperature  distribution  for  a  partly  conducting  crack  ( Bc  =  l) . 


p, 


p, 


- ► 

a 

symmetry  plane 


Figure  2.  Plane  geometry  of  a  thermally  loaded  center  cracked  plate  showing  the  resultant  force  and  moments  at  the  central 
symmetry  plane.  Surface  temperatures  T^T2  are  fixed  with  Tx  >  T2  ,  while  crack  face  temperatures  T^+\t^  are 
dependant  on  position  along  the  crack  face. 
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r/H 


Figure  3.  Piecewise  stress  distribution  function.  Examples  shown  are  for  77  =  1 ,  A  =  1  and  Bc  =  0  .  For  the  symmetric 
case  the  absolute  value  of  the  stress  (shown)  is  the  same  on  upper  and  lower  faces.  The  theory  lines  are  calculated  using 
Eq.  (14)  with  Kn  and  <T0  a)  extracted  from  the  finite  element  solution  (solid  lines)  and  b)  based  in  the  large  a  limit, 

|cr0|  =  Ea{Tx -T2)/2  ,  Kn  =  Ea{Tx—T2^{2Hy2  /$>  (dashed  line). 


a/H 


Figure  4.  Average  temperature  jump  across  crack  face  for  a  centered  crack  (rj  =  l) .  Markers  are  based  on  finite  element 
results,  calculated  using  Eq.  (12).  Theory  lines  are  Eq.  (25)  with  fit  parameter^  =  0.21 12  .  The  fit  is  remarkably  good  for  a 
wide  range  of  Bc  and  a/H  . 
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a/H 


Figure  5.  Energy  release  rate,  for  r/  =  1 ,  A  =  1  and  Bc  =  0  .  Theory  lines  are  Eq.  (21)  using  (j)  =  1 ,  for  a  fully  analytic 
result  that  demonstrates  an  improvement  over  the  infinite  crack  length  solution.  Using  the  improved  value  of  (j)  from  Eq.  (25) 
yields  an  excellent  fit  to  the  numerical  results.  The  calculation  is  remarkably  well  behaved  for  shorter  cracks  than  one  might 
expect  (dashed  region),  it  does  however  become  non-conservative  go  to  zero  at  nonzero  a  . 


Figure  6.  Energy  release  rate  for  the  isotropic  case  (A  =  l)  ,  for  7/  =  1,1/2  ,  and  Bc  =  0,1 .  The  theory,  Eq's.  (21)  and  (25), 
fits  best  for  7j  —  1 ,  and  provides  a  conservative  estimate  for  off-center  cracks.  The  numeric  values  to  the  right  are  the  exact 
asymptotic  values  for  the  four  cases  from  Eqn's  (7),  and  (20). 
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HjEX 


Figure  7.  Orthotropic  example.  Normalized  energy  release  rate  for  a/ H  =  4,/t  =  l/2,2?c  =0  showing  the  variation  with 
shear  modulus.  Reduced  energy  release  rate  is  predicted  by  both  theory  and  numerical  simulation,  but  only  for  very  small 
values  of  /uxz . 
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